Nuclear Structure Calculations with Low-Momentum Potentials in a Model Space 

Truncation Approach 
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We have calculated the ground-state energy of the doubly magic nuclei 4 He, ie O and 40 Ca within 
the framework of the Goldstone expansion starting from various modern nucleon-nucleon poten- 
tials. The short-range repulsion of these potentials has been renormalized by constructing a low- 
momentum potential Vj ow _k- We have studied the connection between the cutoff monemtum A 
and the size of the harmonic oscillator space employed in the calculations. We have found a fast 
convergence of the results with a limited number of oscillator quanta. 
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I. INTRODUCTION 

One of the well-known important features of realistic 
nucleon-nucleon (NN) potentials is their strong repul- 
sive behavior in the high-momentum regime. This im- 
plies that when performing nuclear structure calculations 
within a perturbative approach it is unavoidable to renor- 
malize the NN potential. The renormalization may be 
achieved by constructing a potential that takes into ac- 
count the high momentum components of the original po- 
tential in an effective way. For example, the well-known 
Brueckner reaction matrix G Q is an energy-dependent 
effective interaction, defined in a model space P, obtained 
projecting out all two-particle excitations above a chosen 
Fermi surface. The main drawback of such a procedure 
is that the G matrix is energy dependent, which stems 
from the fact that it does not fulfill a decoupling con- 
dition between the model space P and its complement 
Q. 

Recently, we have proposed 0, a new method to 
renormalize the bare NN interaction, which is proving to 
be an advantageous alternative to the use of the Brueck- 
ner G matrix. We define a low-momentum P-space up 
to a cutoff momentum A and derive from the original 
Vnn an effective low- momentum potential Vjow-k that 
satisfies the decoupling condition between the low- and 
high-momentum spaces. This Viow-k is a smooth poten- 
tial which preserves exactly the on-shell properties of the 
original Vnn and is suitable for being used directly in 
nuclear structure calculations. 

In applying the low-momentum JVJV potential Vj ow _k 
to nuclear structure calculations, an important step is 
the determination of the decimation momentum A for 
VJow-k- There are two considerations: first, A is based 
on the separation of scale idea of the renormalization- 
group effective field theory (RGEFT) approach. In low- 
energy nuclear physics, one is interested in phenomena 
at low-energy scale and consequently the details of the 
short-distance (high-energy scale) physics are unimpor- 



tant. This leads to the derivation of the low- momentum 
potential Vow-k by integrating out the high-momentum 
components of modern NN potentials beyond a cutoff 
momentum A. 

The second consideration is that all modern NN po- 
tentials are constructed to fit the empirical phase shifts 
up to the inelastic threshold Pi aD ~ 350 MeV, which cor- 
responds to a maximum relative momentum ~ 2.1 fm -1 . 

In the past few years, we have profitably employed 
this technique in realistic nuclear structure calculations 
for both doubly closed-shell nuclei, within the framework 
of the Goldstone expansion 0, [j| , and open-shell nuclei 
within the multilevel shell-model framework [SibJIfilHIS]- 
In all these works we have used for the cutoff momentum 
the value A = 2.1 fm _1 . 

It is interesting, however, to investigate the relation 
of the cutoff momentum A to the dimension of the con- 
figuration space in the coordinate representation, where 
actually our calculations are performed. The study of 
such a relation is the aim of the present work, where 
we show how the choice of a cutoff momentum implies 
a maximum value for the energy of the two-nucleon sys- 
tem, the latter introducing a simple criterion to choose 
the two-nucleon model space. 

To verify the reliability of this approach we calculate, 
in the framework of the Goldstone expansion, the ground 
state (g.s.) properties of doubly closed-shell nuclei start- 
ing from different NN potentials renormalized through 
the Viow-k procedure. 

The paper is organized as follows. In Sec. II we give an 
outline of our calculations. A detailed discussion of the 
convergence properties of the perturbative calculations is 
presented in Sec. III. Sec. IV is devoted to the presen- 
tation and discussion of our results for the three nuclei 
4 He, 16 0, and 40 Ca. Some concluding remarks are given 
in Sec. V. 
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II. OUTLINE OF CALCULATIONS 

As already mentioned in the Introduction, we renor- 
malize the short-range repulsion of the bare NN poten- 
tial Vmn by integrating out its high momentum compo- 
nents through the so-called Viow-k approach. A detailed 
description of this procedure can be found in Refs. 
The renormalized NN potential Vj ow -k preserves the ob- 
servables predicted by Vmn f° r the two-nucleon system, 
and consequently its x 2 / datum, up to the cutoff momen- 
tum A. The Viow-k is a smooth potential and can be 
used directly in nuclear structure calculations within a 
perturbative approach. 

While the Vj ow -k is defined in the momentum space, 
we perform our calculations for finite nuclei in the co- 
ordinate space employing a truncated harmonic oscilla- 
tor (HO) basis. Since the Vi ow -k procedure decouples 
the momentum space into the low- and high-momentum 
regime, it is desirable to recover such a decoupling in the 
HO space. 

Let us consider the relative motion of two nucleons in 
a harmonic oscillator well in the momentum representa- 
tion. For a given maximum relative momentum A, the 
corresponding maximum value of the energy is: 



h 2 A 2 
M 



(1) 



where M is the nucleon mass. 

This relation may be re-written in terms of the maxi- 
mum number Af max of HO quanta for the relative coor- 
dinate system, so for a given HO parameter hw we have: 



N, 



h 2 A 2 



(2) 



The above equation provides a simple criterion to map 
out the two-nucleon HO model space. If we write the 
two-nucleon states as the product of HO wave functions 



\ab) = \nJaja, nbhjb) 



(3) 



our HO model space is defined as spanned by those two- 
nucleon states that satisfy the constraint 



this case, we clearly should use this value for A. Let us 
consider another situation where the model-space wave 
functions have important momentum components of, say, 
up to 3.0 fm _1 . Then in this case, we need to use a larger 
A. 
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FIG. 1: Different P n{ 's for 7iw=27 MeV. See text for details. 

The above considerations may be illustrated by in- 
specting the explicit expression of the Vj ow -k matrix el- 
ements in terms of wave functions of the center-of-mass 
and relative coordinates: 

(nJaja, n b l b j b ; JT\V\ OV! -k\nJcjc, n d l d j d ; JT) = (-1) 



L+L' 



LL'S 



.in.li.j, .hi^'-L + I ){'2L' + I ) ! 2 .S" -i- II { 

L S J 

[l-{-l) l+s+T }(nlN£, L\n a l a n b l b , L) 

nln'VNC 

(n'l'NC, L'\nJ c n d l d , L') ^(2J + 1) [ j ' j ^ 
(nlSTJ\V low ^\n'l'STJ) , 




C I' V 
S J J 



(5) 



where j = + 1, (nlNC, L\n a l a n b l b , L) are the 
Brody-Moshinsky transformation brackets [To| . and, ac- 
cording to (4), 2n + 1 and In! + 1' are both < N max . The 
matrix element {nlSTJ\V\ OV! -k\n'VSTJ) is expressed in 
terms of the momentum-space HO wave functions P n is 



as 



+ 2n b + l b < A max . (4) 

The need for such a boundary condition for our model 
space may be also pointed out by the following consid- 
erations. The momentum contents of the two-nucleon 
wave function clearly depend on the choice of the model 
space. For example, the average momentum of HO wave 
functions is proportional to (hu) 1 ^ 2 . Suppose we are do- 
ing a calculation using a small model space together with 
a small fku, so that the momentum contents of the ba- 
sis wave functions are practically all below 1.5 fm -1 . In 



(nlSTJ\Vi ow - k \n'l'STJ) = / dk dk' kk 1 P nl {k) 

Jo Jo 

Pn>l>{k')V<i S X J {k>k>) . (6) 

Because in Eq. (6) the integrals are evaluated up to 
A, it is desirable to throw away those P n i's which ex- 
tend significantly above the cutoff momentum, and per- 
tain therefore to the high-momentum regime. We have 
verified that applying the constraint (4) amounts to ne- 
glect all P„z's which extend their tail more than ~ 5% 
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above A. Fig. 1 is an explanatory picture where we plot, 
for a given tuo =27 MeV, three momentum space HO 
wave functions P n i's with 2n + I =5 , 7, and 9, respec- 
tively, as functions of k. The vertical dashed line denotes 
a value of A = 2.35 fm _1 , corresponding to N max =7. 

Relation (4) has a general validity, it should be applied 
every time the Vj ow _k matrix elements arc calculated in 
the HO basis. 

In this paper, making use of the present approach, we 
have studied the ground state properties of doubly closed- 
shell nuclei within the framework of the Goldstone ex- 
pansion More explicitly, starting from the purely 
intrinsic hamiltonian 



Pi Pj 

MA 



(7) 



where Vij stands for the renormalized Vnn potential 
plus the Coulomb force, we construct the Hartree-Fock 
(HF) basis expanding the HF single particle (SP) states 
in terms of HO wave functions. The HF basis is then 
employed to sum up the Goldstone expansion, including 
contributions up to fourth-order in the two-body inter- 
action. 

Our calculations are made in a truncated model space, 
whose size is related to the values of the cutoff momen- 
tum A and the fku parameter. The calculations are per- 
formed increasing the A^ max value (and consequently A) 
and varying the ftw value until the dependence on N max 
(A) is minimized. 
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FIG. 2: He ground state energy as a function of the number 
of harmonic-oscillator major shells. The employed two-body 
interaction is a VW-k derived from CD-Bonn potential with 
A = 2.1 fm . The oscillator parameter Kuj is equal to 18 
MeV. Calculations have been performed with different Pade 
approximants, the red line representing an exact calculation. 



III. CONVERGENCE PROPERTIES OF THE 
PERTURB ATIVE EXPANSION 

In section IV we will show the calculated ground state 
energies of some doubly closed-shell nuclei within the 



framework of the perturbative approach (Goldstone ex- 
pansion). Here, we study the convergence properties of 
the perturbative series. To this end, we have performed 
a test calculation, starting from a Vj ow -k with a fixed 
cutoff momentum A = 2.1 fuC 1 and derived from the 
NN CD-Bonn potential 12]. For this Vi ow _k, hermitized 
according to the procedure based on the Cholesky de- 
composition suggested in Ref. |13| . an exact calculation 
of the ground state energy of 4 He has been performed 
in the framework of the hyperspherical harmonic (HH) 
approach The value obtained is -27.95 MeV 

considering the Viow-k as a NN potential defined in the 
whole momentum space, with View— k(k,k') = when k 
or k' > A. With the above potential we have calculated 
the same quantity in the framework of the Goldstone ex- 
pansion. 

Using Pade approximants |l6l Il7| one may obtain a 
value to which the perturbation series should converge. 
We consider the following three Pade approximants: 



[L\1] = E + E 1 
where L = 2, 3, and 



+ 



E L 



1 — El+i/El 



[2|2] 
where 

7i = 



E (l + 7i + 72) + £q(l + 72) + E 2 
I + 71+72 



(8) 



(9) 



E2E4 — E 3 
E1E3 



El 



72 = 



E3 + £1 71 
E~ 2 



Ei being the zth order energy contribution in the Gold- 
stone expansion. 

In Fig. 2 the 4 He calculated ground state energy is 
plotted versus the number of HO major shells included in 
the calculation, the red line representing the exact value. 
It is worth here to mention that no further truncation 
for the two-nucleon states has been performed, we con- 
sidering the Vjow-k as defined in the whole momentum 
space. 
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FIG. 3: 4 He ground state energy with CD-Bonn potential as 
function of iV max , for different values of Hu>. The straight line 
represents the Faddeev-Yakubovsky result. 
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FIG. 4: Same as Fig.3, but for N 3 LO potential. 

We see that the Pade approximants [3|1] and [2|2] give, 
for a given number of HO major shells, almost the same 
value , the difference being at most 45 keV. The Pade ap- 
proximant [2|1] differs at most by 200 keV from the [3|1], 
and 245 keV from the [2 1 2] . Moreover, the results, cor- 
responding to the largest space we have employed, come 
close to the exact value, the energies being -27.79 and 
-27.84 MeV with the [3|1] and [2|2] approximants, re- 
spectively. 

On these grounds, we report in the following section 
the results obtained using the Pade approximant [2 1 2] . 
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IV. RESULTS 

To test the reliability of our calculations we have cal- 
culated the binding energy of 4 He starting from different 
Vjvjv's, and compared our results with those obtained 
through the Faddeev-Yakubovsky (FY) method. 

In Figs. 3, 4, and 5 we show the calculated 4 He ground 
state energies obtained from the CD-Bonn N 3 LO 
[l8j. and Bonn A 0] NN potentials, respectively. In 
eac h fig ure the straight red line indicates the FY result 
[20l l2l| while the other curves represent our calculated 
values, for different values of fku, versus the maximum 
number of HO quanta -/V max that binds the two-nucleon 
configurations according to the relation (4). In Fig. 6, 
we report the same results of Fig. 3, but versus the cutoff 
momentum A. 

For the CD-Bonn, N 3 LO, and Bonn A potentials we 
obtain convergence with hw — 36, 29, and 34.5 MeV, 
respectively. Our calculated energies are -25.92, -25.02, 
and -27.78 MeV for the above Vm n 's. These values are in 
good agreement with the FY results, the largest discrep- 
ancy being 0.39 MeV for N 3 LO potential. We have done 
similar calculations starting from other modern NN po- 
tentials, such as the Nijmegen II [22 and Argonne V18 
[23j potentials, but because of their stronger tensor com- 
ponents it has not been possible to achieve convergence. 
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FIG. 5: Same as Fig.3, but for Bonn A potential. 



FIG. 7: O ground state energy with CD-Bonn potential as 
function of Af max , for different values of Hu>. The straight line 
represents the experimental value |25l| . 
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FIG. 6: 4 He ground state energy with CD-Bonn potential as 
function of cutoff momentum A, for different values of hui. 
The straight line represents the Faddeev-Yakubovsky result. 



We have calculated also the ground state energies of 
16 and 40 Ca starting from the CD-Bonn potential. In 
Fig. 7 we present our results for 16 O and compare them 
with the experimental datum. The converged value, as 
can be seen from the figure, is obtained for Hui = 27 
MeV and is equal to -118.1 MeV, the discrepancy with 
the experimental value being 9.5 MeV. In this case we 
cannot compare our calculations with the exact ones. 
It is worth mentioning, however, the work by Fujii et 
al. [24| who, using the unitary model-operator approach 
(UMOA), predict for the CD-Bonn potential a converged 
value of -115.61 MeV, including only two-body correla- 
tions. In a more recent paper [26|. the above authors 
have estimated the three-body cluster effect to contribute 
about -4 MeV. 
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As regards 40 Ca, a calculation of its ground state en- 
ergy including fourth-order contributions in the Gold- 
stone expansion has not been done because of the large 
CPU time needed. We therefore report in Fig. 8 the 
results obtained with the Pade approximant [2|1], taking 
into account up to third-order contributions in the Gold- 
stone expansion. The converged value is -307.8 MeV with 
huj = 25.5 MeV. In this case, the discrepancy with respect 
to the experimental value is 34.2 MeV. 



V. SUMMARY 

In this work, we have calculated the ground state en- 
ergy of some doubly closed-shell nuclei in the frame- 
work of the Goldstone expansion, starting from differ- 
ent realistic NN potentials. The short-range repulsion 
of these potentials has been renormalized by integrating 
out their high-momentum components through the so- 
called Viow-k approach. We have introduced a criterion 
to map out the model space made up by the two-nucleon 
states in the HO basis, according to the value of the cutoff 
momentum A. 



-280 

„ -300 
> 

CD 

2, -320 

LU 

-340 
-360 



CD-Bonn Ca 



• ft co=26.5 

• 7[q3=26_ 

• H co=25.5 

• H 03=25 

• n 03=24.5 




34.2 MeV 



8 

N„„ 



10 



12 



14 



FIG. 8: Same as Fig.7, but for for 40 Ca 



The reliability of this procedure has been tested by 
calculating the ground state energy of 4 He, with the CD- 
Bonn, N 3 LO, and Bonn A potentials and comparing the 
results with the FY ones. We have found that the energy 
differences are at most 390 keV. These differences are 
due to two reasons. On the one hand, our calculations 
have been performed using a perturbative approach, so 
that small contributions coming from higher order terms 
may have not been completely taken into account by the 
Pade approximants. On the other hand, we do not expect 
that Eq. (2) recovers exactly in the HO basis the V[ w-k 
decoupling into low- and high-momentum regime. 

In any case, the limited size of the discrepancies shows 
that our approach provides a reliable way to renormalizc 
the NN potentials preserving not only the two-body but 
also the many-body physics. 

On the above grounds, we have performed similar cal- 
culations for heavier systems, such as 16 and 40 Ca and 
obtained converged results using model spaces not ex- 
ceeding N max = 10. 

The rapid convergence of the results with the size of 
the HO model space makes it very interesting to study 
in a near future heavier systems employing our present 
approach. 
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